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We present a detailed study of the effect of local chemical ordering on the structural, electronic, 
and dynamical properties of amorphous gallium arsenide. Using the recently-proposed "activation- 
relaxation technique" and empirical potentials, we have constructed two 216-atom tetrahedral con- 
tinuous random networks with different topological properties, which were further relaxed using 
tight-binding molecular dynamics. The first network corresponds to the traditional, amorphous, 
Polk-type, network, randomly decorated with Ga and As atoms. The second is an amorphous struc- 
ture with a minimum of wrong (homopolar) bonds, and therefore a minimum of odd-membered 
atomic rings, and thus corresponds to the Connell-Temkin model. By comparing the structural, 
electronic, and dynamical properties of these two models, we show that the Connell-Temkin net- 
work is energetically favored over Polk, but that most properties are little affected by the differences 
in topology. We conclude that most indirect experimental evidence for the presence (or absence) of 
wrong bonds is much weaker than previously believed and that only direct structural measurements, 
i.e., of such quantities as partial radial distribution functions, can provide quantitative information 
on these defects in a-GaAs. 
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I. INTRODUCTION 

After twenty five years of effort, the structure of amor- 
phous materials, and how it affects the electronic and 
vibrational properties, remains largely unresolved. Most 
experimental probes yield information that is averaged 
out over rather large lengthscales and therefore lack the 
sensitivity to discriminate between various possible struc- 
tural models of the same material. Techniques such as 
EXAFS, while they can provide structural information 
at the atomic level, are often too imprecise to yield defi- 
nite and unambiguous structural parameters. One must 
therefore proceed iteratively between models and exper- 
imental data in order to acquire the desired structural 
information. 

In the case of amorphous semiconductors, progress in 
the development of satisfactory structural models has 
been hindered by difficulties in constructing "continuous 
random networks" (CRN) with different topologies, i.e., 
appropriate to different materials. The idea of represent- 
ing the structure of amorphous semiconductors by CRN's 
was first proposed by Zachariasenp in this picture, the 
material is assumed to consist of a "collage" of tetrahedra 
quite similar to those found in the corresponding crystal 
but randomly connected through their vertices. Based 
on these ideas, the mechanical CRN constructed by Polk 
was found to provide a very satisfactory descriptionj-pf 
the topology of elemental amorphous semiconductorsJa 

For compound materials, however, the situation is not 
as clear: The building-block tetrahedra, if they exist. 



can be formed in many different ways, depending on 
the chemical identities of the atoms, the arrangement of 
which is determined by the bonding characteristics of the 
material. For example, in the case of the IIPV compound 
GaAs, chemical ordering should predominate because the 
material is partly ionic, i.e., the number of bonds between 
like atoms ( "wrong bonds" ) should be minimal. Ideally, 
a tetrahedron should consist of a Ga atom surrounded by 
exactly four As atoms (or vice versa), as is the case in a 
perfect (zinc-blende) GaAs crystal. The actual structure 
of the disordered material, therefore, will be determined 
by a balance between the cost of elastically deforming 
the network while maintaining perfect chemical ordering 
and the cost of introducing wrong bonds. In view of this, 
a-GaAs appears to be an ideal camlidate for the realiza- 
tion of the Connell-Temkin modefl — a CRN similar to 
Polk's but without odd-membered atomic rings. (A ring 
is defined as a closed path between an atom and itself 
through a series of bonds). Indeed, in an unconstrained 
CRN, there are inevitably both even- and odd-membered 
rings. Since odd-membered rings necessarily bring about 
wrong bonds, which cost Coulomb energy, it is expected 
that the number of them will be minimal in a-GaAs, and 
ideally none as in the Connell-Temkin model if the cost 
in elastic deformation energy associated with this con- 
straint is not too large. 

As noted above, detailed experimental information re- 
garding the structure of compound materials is in gen- 
eral difficult to obtain; this is particularly so in the case 
of a-GaAs because of the close similarity between the 



constituent elements. (They are near neighbors in the 
Periodic Table; in fact, the similarity in size is another 
reason why the topology of a-GaAs should correspond 
closely to that of the single-component Connell-Temkin 
model.) Experimental evidence for the presence (or ab- 
sence) of wrong bonds in this material is indeed essen- 
tially non-existent, while the question does not arise in 
elemental amorphous semiconductors, such as a-Si. It is 
therefore important, in order to understand on a local 
scale the topology of amorphous semiconductors, to ex- 
amine idealized representations of the two different types 
of networks, viz. Connell-Temkin and Polk, and compare 
them with experiment. This is the object of the present 
paper. 

In a recent article, we have shown that it is quite 
possible to construct a binary-compound network which 
contains essentially no odd-membered rings (i.e., wrong 
bonds)p. corresponding to an infinite Connell-Temkin 
model;Q for GaAs, this model is found to be energeti- 
cally preferred over the Polk model, where both even- 
and odd-membered rings are present. Likewise, for ele- 
mental and/or non-ionic tetrahedral semiconductors, the 
Polk- type model is preferred for entropic reasons: elas- 
tic costs associated with the constraint on the absence 
of odd rings are negligible but the limitation on rings 
restricts significantly the number of possible configura- 
tions. In this way, we have been able to achieve a direct 
structural comparison of materials differing on the inter- 
mediate lengthscale — a-Si and a-GaAs. 

In the present paper, we extend this study and exam- 
ine in detail the structural, vibrational, and electronic 
properties of amorphous GaAs as described by the two 
types of CRN's mentioned above. This provides unique 
and much-needed information on the effects of topology 
on the properties of amorphous semiconductors. Indeed, 
upon comparing different networks constructed using the 
same energy scheme, it is possible to isolate, in measured 
properties, those effects arising from topology from those 
due to the interactions between atoms — evidently some- 
thing which cannot be done experimentally. 



II. METHODOLOGY 

The timescale on which chemical ordering takes place 
when a binary compound is cooled down from the liquid 
phase depends strongly on the ionicity of the material. 
For example, molecular-dynamics (MD) simulations on 
a timescale of picoseconds are sufficient to ensure proper 
ordering of silica (Si02) p which has a Phillips ionicity of 
2.09.B In contrast, corresponding simulations of a-GaAs, 
which is only slightly ionic (0.22), have not given clear 
indication of chemical ordering. Whether these results 
reflect the actual structure of GaAs or limitations of MD 
simulations needs to be addressed using a different ap- 
proach for constructing structural models. 

One possibility is to bypass the dynamics of formation 



and devise an appropriate static structure optimization 
scheme; by "appropriate" we mean which will lead to a 
physically realistic structural model. The route followed 
need not be physical, however: the philosophy we adopt 
here is one where "the end justifies the means". This 
argument will become evident in the discussion which 
follows, but one closely-related precedent can be invoked 

— the celebrated Wooten- Winer- Weaire (WWW) algo- 
rithm for constructing models of a-Si:Ll In this approach, 
crystalline Si is amorphized through a sequence of bond- 
switching moves which are totally unphysical; yet, the 
final (converged) structure possesses much of the prop- 
erties of real a-Si. The method cannot be employed in 
the case of compound semiconductors because it is the 
essence of it to introduce odd-membered rings and wrong 
bonds. (It involves breaking nearest-neighbor bonds and 
forming second- neighbor — thus wrong — bonds.) These 
limitations of the WWW algorithm, and/or absence of al- 
ternative models, have significantly hindered the study of 
amorphous semiconductors over the last ten or so years. 

The activation-relaxation technique (ART) recently 
proposed by Barkema and Mousseauj providesr-a way to 
circumvent the restrictions of the WWW modelB Given a 
model of interatomic potentials, which can be of any form 

— from purely empirical to fully ab initio, the method 
(which is event-based) forces relaxation through a series 
of physical moves in a configurational space reduced to 
a set of isolated energy minima connected together by 
paths going through saddle points on the configurational 
energy landscape. Since the moves are defined in the 3A^- 
dimensional configurational space, where N is the num- 
ber of atoms, the algorithm is completely independent of 
the structure in the 3-diniensional real space. A move, 
therefore, can involve any number of atoms, and in par- 
ticular be local or span the whole system; it is not limited 
to a pre-determined list of events, such as bond switches 
d la WWW or atomic exchanges. Further, since the con- 
figurational space is reduced to an (infinite) number of 
discrete points, events can be defined uniquely. 

Full details and illustrations of the ART method can 
be found in Ref. j^; here we give a brief overview. An ART 
simulation starts with the system in a local minimum of 
the potential-energy surface. The configuration is then 
placed slightly out of equilibrium by operating a small 
change on the system, e.g., moving an atom in a random 
direction by a very small amount. The component of the 
force parallel to the displacement of the configuration is 
then inverted and the whole configuration pushed away 
from the nearby minimum following the modified force: 



G = F - [1 + a)iF ■ AX)AX, 



(1) 



where G and F are both 3A''-dimensional vectors, a is a 
dimensionless parameter (here set to 0.15) and X is the 
unit-vector displacement from the local minimum to the 
current position. Equation ms iterated until the modified 
force (and thus also the real force) vanishes, indicating 
that a saddle point has been reached. The configuration 



is then pushed sHghtly away from the saddle point and 
relaxed to a new local minimum. This procedure is iter- 
ated until convergence, i.e., until no further changes in 
the total energy are observed. 

A few remarks are in order: (i) Since the algorithm is 
defined in configurational space, there is no limit to the 
number of atoms which can be involved in a single event; 
for the materials studied here, events involving one to a 
few tens of atoms have been observed, (ii) In order to 
prevent the configuration from moving to close-by saddle 
points which are not significant at finite temperature, a 
repulsive force is introduced that excludes a region of 
width r^ Xc around any local minimum: 



^rep — r^\X ^cji 



(2) 



here x is the scalar displacement of the configuration from 
the local minimum, Xc is a cutoff parameter, and A is the 
strength of the repulsive part. Both x^. and A are drawn, 
for each new event, from a linear random distribution: 
0.3 <Xc< 1.3 A and 0.5 < A < 1.5 eV/A^. (in) After 
each event, the volume of the system is optimized such 
as to minimize the configurational energy, (iv) Both ac- 
tivation to the saddle point and relaxation to the near- 
est local minimum are-jDcrformcd using the Lcvenberg- 
Marquardt algorithm;ty this algorithm, which includes 
both the steepest-descent and the Hessian approaches, 
is fairly efficient around a minimum and does not misbe- 
have far from it. (v) The optimization is performed using 
a multiple-configuration simulated-annealing approach; 
thus, each new event is accepted with a probability 



Paccept{{xi+l}) = exp 



E({x,+i}) - E{{x,}) 



keT 



, (3) 



where here the temperature is a is non-physical param- 
eter. In practice, two configurations are run in parallel 
at different temperatures. After each step, the energies 
are compared and the configurations are switched accord- 
ing to a Metropolis rule. The use of multiple configura- 
tions allows for more efficient sampling of configurational 
space, permitting configurations to escape from dead-end 
minima, i.e., those which cannot lead directly to a lower- 
energy state. 

The force F in Eq. n^ is derived from an interatomic po- 
tential which, as noted above, can in principle be of any 
form. Evidently, the final, optimized structure will de- 
pend on the choice of potential. With respect to GaAs, 
there exists (to our knowledge) no satisfactory empir- 
ical potential for the disordered phase. Our own at- 
tempts in thisojegard using a recently-proposed Tersoff- 
type potentialtll have lead to structures in deep disagree- 
ment with experiments. For the present work, construc- 
tion of the computer models proceeded in two stages un- 
der two sets of potentials: first, ART optimizations us- 
ing modified Stillinger- Weber potentials (see below) were 
carried out; and second, the models were further relaxed 
under semi-empirical tight-binding (TB) potentials. The 
reason for this "double-relaxation" approach is that while 



ART relaxations can in principle be done with the TB 
potentials, this remains a complicated and computer- 
intensive enterprise. Carrying out a "first pass" with 
empirical potentials allows optimized structures to be ob- 
tained rapidly, while the final TB relaxation provides a 
physically-meaningful basis for the models; indeed, we 
have found the properties of our networks to be only lit- 
tle affected by the final TB relaxation, thus indicating 
the convergence and validity of the procedure. 



III. MODEL PREPARATION 

As discussed in the Introduction, we consider here 
two networks with different topologies; following Ref. |j, 
these are poetically labeled CRN- A and CRN-B. CRN- 
A corresponds to a Polk-type network, while CRN-B is 
Connell-Temkin like. The initial ART relaxation was per- 
formed, in_both cases, with a Stillinger- Weber potential 
for silicon Jl3 except that, in order to compensate for toi 
weak an angular force in the original set of parameters,! 
the three-body contribution was increased by 50%. For 
CRN-B, moreover, a repulsive term between like atoms 
was introduced in order to minimize the number of wrong 
bonds and favor a Connell-Temkin like topology: 
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(4) 



where e is the energy parameter of the Stillinger- Weber 
potential, A^j = 1.2 for like atoms and zero otherwise, 
and s = 3.6 A; this is also the value for the cut-off of this 
potential. Except for this additional repulsive energy, 
both atomic species are treated in exactly the same way 
for CRN-B. The validity of these modifications is assessed 
a posteriori when the lattice is further relaxed with a 
physical TB potential (see below). As with any static 
approach to modeling complex systems, it is the final 
structure which determines the quality of the method, 
i.e., here, the end justifies the means, as noted earlier. 

The ART optimization was initiated, in each case, from 
a 216-atom random-packed configuration at the c-Si den- 
sity; periodic boundary conditions were employed in or- 
der to eliminate surface effects. Iteration of the ART pro- 
cedure was carried out until changes in the total energy 
were judged insignificant. For CRN-A, the optimization 
is fairly rapid and takes two or three days on a R8000 
computer. In the case of CRN-B, considerable atomic 
diffusion is required in order to minimize the number of 
wrong bonds, and the simulation takes longer — about 
one week. These run times should be contrasted with 
the several months needed for a TB simulation p£.a 64- 
atom system on a workstation or a Car-Parrinelloll^ (CP) 
simulation on a state-of-the-art parallel supercomputer. 

Both CRN-A and CRN-B were then relaxed with both 
the, Goodwin-Skinner-Pcttifor (GSP) TB potential for 
SillJ and the Melteni-Colombo-Migfio (MCM) TB poten- 
tial for GaAs.t^ Thus, in total, we have four different 



zero-temperature TB-relaxed models: CRN-A-Si, CRN- 
A-GaAs, CRN-B-Si, and CRN-B-GaAs, i.e., for each ma- 
terial, two models with different topologies. In the case 
of CRN-A-GaAs, it is necessary to "label" the atoms af- 
ter the ART simulation and before the TB relaxation. 
This was done using a random "label-switching" proce- 
dure designed to minimize the number of wrong bonds 
on the lattice. In this way, we obtained a proportion of 
wrong bonds of 14%, as close as one can possibly get, for 
a finite-size system, to the theoreticalr-value of 12% for 
optimal ordering on a Polk-type CRN.E3 The configura- 
tional energies of each of the four models after the static 
TB relaxation are presented in Table |. A detaiied com- 
parison between Si and GaAs is given elsewhereQ and we 
therefore only discuss, here, the two models for GaAs, 
viz. CRN-A-GaAs and CRN-B-GaAs, hereafter referred 
to simply as CRN- A and CRN-B. 

Following the static TB simulation, the two samples 
were further relaxed at 300 K using MD for a total of 
7.0 ps. Although this is a fairly short period of time, 
it is enough to ensure that both CRN-A and CRN-B 
have evolved into deep local minima of the potential en- 
ergy surface. In order to verify this, we have also an- 
nealed the samples at 700 K during 8.8 ps before run- 
ning again at 300 K for 3.5 ps and 10 K for 0.9 ps. The 
CRN-B network was found to be only very weakly af- 
fected by the high-temperature annealing. For CRN-A, 
in contrast, annealing resulted in significant changes in 
the topology; the average coordination, for instance, in- 
creased from 3.95 to 4.19. Likewise, the energy increased 
from —13.45 eV/atom to —13.39 eV/atom. This in fact 
could be expected: it indicates that model CRN-A is not 
a proper state for GaAs. Though it would take much 
longer runs to find out, it is very likely that the system is 
trying to find a route (through a higher-energy transition 
state) towards the preferred configuration — one without 
wrong bonds. Thus, these high-temperature simulations 
confirm that CRN-B is indeed a much better model for 
the structure of a-GaAs than CRN-A. In the discussion 
that follows, the structures relaxed at 300 K before an- 
nealing will be used for comparing the two models. 

As noted above, here we used 216-atom unit cells, the 
largest size we could deal with in our TB simulations, 

TABLE I. Energies in eV per atom of the two networks 
statically relaxed using the two sets of tight-binding param- 
eters (see text). The atoms on the CRN-A model have been 
placed in such a way as to minimize the number of wrong 
bonds. Also given are the energies of the crystalline phases. 
SL refers to the tight-binding simulation of Seong and Lewis, 
Ref. I g^ 
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FIG. 1. Ball-and-stick representation of CRN-A and 
CRN-B after relaxation with a GaAs TB potential; small 
white circles are As and larger grey circles are Ga atoms. 

nevertheless large enough to provide a satisfactory de- 
scription of amorphous semiconductors. With 64- atom 
cells, we found ART to lead, in all cases, to the crys- 
talline state, indicating that the configurational space for 
a system of this size is small enough for ART to find the 
global minimum. In contrast, for the 216-atom systems, 
crystallization never occured and so we conclude that the 
local minimum we find are truly optimized. 



IV. STRUCTURAL PROPERTIES 

In Fig. |l| are shown ball-and-stick representations of 
the two a-GaAs samples. For CRN-B, wrong bonds 
are few but one can be seen in the top right quadrant 
where its presence gives rise to both a five- and a seven- 



membered ring. It is clear from this figure that the two 
models, though clearly ordered at short-range — both 
topologically and chemically — bear no trace of crys- 
tallinity. (It is often the case that partly crystalline sam- 
ples cannot be recognized, because of averaging and ther- 
mal agitation, in such quantities as radial distribution 
functions and static structure factors.) In fact, based on 
this visual inspection, it is already evident to a trained 
eye that the samples are excellent models of the amor- 
phous material. 

We now proceed with a more quantitative analysis 
of the models. First we give, in Table ||, some of the 
usual system-averaged structural quantities — coordina- 
tion numbers, proportion of wrong bonds, and width of 
the bond-angle distribution — for our two models, both 
at and at 300 K. Here we do not distinguish between 
the various types of correlations and treat all atoms on 
the same footing. The reason for this is that the identity 
of the atoms prior to the ART relaxation was introduced 
in an ad-hoc fashion; since the TB relaxation does not 
change the connectivity of the network, the two species 
are still topologically equivalent in the final configura- 
tion. This, it turns out^ consistent with previous simu- 
lations by Molteni et altl and by Fois et alE3 which show 
the two species to behave symmetrically, and is also sup- 
ported by experiment, which shows As and Ga to be both 
four- fold coordinated (taking due account of variations in 
composition) .E3 

We first observe that the structural characteristics of 
both CRN-A and CRN-B satisfy the requirements of a 
"good" CRN, namely four-fold coordination and small 
bond-angle deviation. From Table ||, it is apparent that 
CRN-A and CRN-B have a density of coordination de- 
fects lower than that of previous models of a-GaAs. In 
particular, both models have a coordination of almost 
exactly four, i.e., most atoms are perfectly coordinated 
but a few are under-coordinated (cf. Table ||) . [In order 

TABLE II. Structural characteristics of the models dis- 
cussed in the text, at both K and 300 K: distribution of co- 
ordination numbers, Z (and nearest-number cutoff distance, 
rjvjv), density of wrong bonds, and width of the bond-angle 
distribution, AO. Also given are the corresponding numbers 
from other simulations, all at K: SL — tight-binding simu- 
lations of Ref. Ell; MCM — tight-binding simulations of Ref. 
|l7|; CP — Car-Parrinello simulations of Ref. [18. 
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FIG. 2. Total static structure factors for CRN-A (dashed 
line) and CRN-B (solid line); the dots are the experimental 
data of Udron et ai, Ref. E2l 

to define nearest neighbours, we use here the distance 
corresponding to the minimum following the first maxi- 
mum of the total radial distribution function (see below), 
viz. 3.0 A.] Both networks, therefore, are consistent with 
experiments.E3 It is also clear from Table O that the 
finite-temperature models are totally equivalent to the 
zero-temperture models: Thermal agitation only brings 
about a widening of the distributions of neighbours and 
bond angles, leaving essentially unchanged the number 
of wrong bonds. 

The density of wrong bonds in model CRN-A is 14 or 
15%. As discussed earlier, this value was obtained by 
assigning the identities of atoms on CRN-A such as to 
minimize the number of wrong bonds, and corresponds 
quite closelii,to the "theoretical limit" of 12% for a Polk- 
type CRN.tj It is also close to the values obtained in 
melt-and-quench MD simulations of a-GaAs (10-13% — 
cf. Table H). In contrast, CRN-B, with less than 4 or 5% 
of wrong bonds, is the closest realization of an "infinite" 
Connell-Temkin model. Such a low proportion of wrong 
bonds also seems to be in much better agreement with 
experiment, which indicates |tiiat at most a few percent 
of wrong bonds are present;c3 as will be discussed be- 
low, however, the measurements on which this estimate 
is based turn out to be much less sensitive to the density 
of wrong bonds than what is usually believed. Finally, 
visual inspection of Fig. n reveals no spatial concentra- 
tion of wrong bonds; they appear to be homogeneously 
distributed on both networks. 

That CRN-B is a better model for the structure of 
a-GaAs than CRN-A is evident from the configurational 
energies given in Table ||; it is also evident from this Table 
that the procedure employed here — ART plus TB-MD 
relaxation — leads, in this case, to a better model than 
the usual melt-and-quench MD approach. (Cf. the TB- 
MD results of Seong and Lewis, Ref. El^, given in Table 

I) 

The ability of CRN-B to describe a-GaAs can also be 

assessed from the total static structure factors (SSF's), 
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FIG. 3. Partial static structure factors for the two models, 
as indicated; the dotted, dashed, and solid lines are for the 
Ga-Ga, Ga-As, and As-As partial correlations. 

presented in Fig. 0: the SSF far CRN-B matches more 
closely the experimental dataE3 than that for CRN-A. 
Unfortunately, the samples used in the experiment suffer 
from some inhomogeneities; further structural measure- 
ments on better-quality material would provide much- 
needed experimental data for more accurate comparisons. 

The partial SSF's for the two models are shown in Fig. 
0. Although they differ in many ways, no evident sig- 
nature of the presence of wrong bonds in CRN-A can 
be identified: the differences between the two sets of 
curves are essentially quantitative, and no peaks appear 
in one of the partials which is totally absent in the other. 
However, as we discuss next, the partial radial distribu- 
tion functions (RDF's), obtained from the SSF's by a 
Fourier transformation, allow a much better interpreta- 
tion of these differences. 

In effect, the partial RDF's can provide direct, quan- 
titative, evidence of the existence, and proportion, of 
wrong bonds; they are given in Fig. for our models. 
The fact that CRN-B is chemically ordered is clearly vis- 
ible in the partial RDF's: The unlike-atom partial func- 
tion, Ga-As, exhibits a strong first-neighbour peak, but 
very little amplitude at the second and fourth nearest- 
neighbour distances. In contrast, the like-atom partial 
RDF's, Ga-Ga and As-As, have essentially no nearest- 
neighbours, but exhibit strong second and fourth nearest- 
neighbour peaks. Thus, chemical-order "filters out" the 
shell structure of the material (on the short and inter- 
mediate lengthscales). As a result, the large split peak 
of the Ga-As correlation function in the range 3.5-7 A, 
which corresponds to third ncarcst-ncigbbours, can be 
clearly isolated from others. This is important because 




FIG. 4. Partial radial distribution functions for the two 
models, as indicated; the dotted, dashed, and solid lines are 
for the Ga-Ga, Ga-As, and As-As partial correlations, re- 
spectively. The lower panel gives the total (unweighted) radial 
distribution function for CRN-A (dashed line) and CRN-B 
(solid line). 

this peak corresponds to the various possible dihedral 
conformations; here we find that the two sub-peaks cor- 
respond to dihedral angles of 60 and 180 degrees, as we 
indeed find below through a direct calculation of the di- 
hedral angle distribution. 

These results points to the importance of measuring 
the partial RDF's (or SSF's). Because Ga and As are 
close to one another, however, only the total RDF is 
available from x-ray or othes scattering measurements. 
This is true also of EXAFSl23 for which it is difficult to 
distinguish the atomic type in the nearest-neighbor shell. 
This is unfortunate because, as can be seen in Fig. 0, 
most of the information on wrong bonds is lost in the 
weighted sum over the partial RDF's: the total RDF's 
for the two models are almost identical over the whole 
range of distances of interest, as can be seen in Fig. H. 

The quality of the models can also be inferred from the 
distributions of bond and dihedral angles. In the case of 
a-Si, the width of the bond-angle distribution at K i£ 
found experimentally to lie in the range 10-12 degrees^ 
in accord with recent, fully-optimized, WWW models.B^ 
From MD simulations of a-Si, and more recently a-GaAs, 
typical values for this quantity are in the range 15-17 de- 
grees; in the case of a-GaAs, further, the distobution of 
bond angles is observed to be asymmetric JlJ'Ell biased 
towards smaller angles, manifest of the presence of a sig- 
nificant number of over-coordinated atoms. These results 
could be seen as supporting the analysis of Connell and 
Temkin who have found their model to have a wider 



c 

3 
U 

vS 

c 
o 




60 90 120 150 



3 


(deg.) 


■b> 




IQ 




•o 






^^^U /^ 


J3 


T /\ / - 


(0 


\ J V V " 


J3 


-V^ Nil — ^ - 


O 




u 





QU 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




60 120 180 




(deg.) 





15 




10 


xn 

3 


5 







■« 


10 


Q 





-n — I I I I I — I I I I I — I I I I I I I 






-15 -10 -5 

E (eV) 



FIG. 5. Distributions of (a) bond and (b) dihedral angles 
for CRN- A (dashed lines) and CRN-B (full lines). 

distribution of bond angles than Polk, owing to the ad- 
ditional constraint on the parity of atomic rings. How- 
ever, we find here that CRN-A and CRN-B both have 
a bond-angle-distribution of width about 11 degrees (at 
K; cf. Table O), very much in agreement with experi- 
ment. Both distributions, moreover, are centered closely 
on the ideal tetrahedral angle and are almost Gaussian. 
The similarity between the angular distributions for two 
such different models reflects the ability of the network 
to reorganize in spite of topological constraints on the 
formation of odd-membered rings .c3 

We now turn to dihedral angles, i.e., angles between 
second-nearest-neighbor bonds. Connell and Temkin 
found, in their model, staggered configurations (0 = 60 
degrees) to be four times more numerous than in the Polk 
model, concluding that this preference for staggered con- 
figurations should be a signature of the absence of odd- 
membered rings, i.e., characteristic of CRN-B-type ma- 
terials. We find no support for such a conclusion here: 
as can be seen in Fig. a(b), our two models exhibit es- 
sentially identical dihedral-angle distributions. 

As demonstrated in Ref. 0, a-GaAs and a-Si form net- 
works which are topologically different. Structural signa- 
tures of these differences, however, as we have seen here 
by comparing models GRN-A and GRN-B, are extremely 
difhcult to extract from experiment or even from com- 
puter models; it seems to show up clearly, in fact, only 
in quantities which cannot be measured directly, namely 
the number of wrong bonds or ring statistics. Thus, 
most measurable quantities appear to be unaffected by 
the constraint on wrong bonds. It should be noted that 
even though the experimental precision required to de- 
cide between the two models on the basis of their total 



FIG. 6. Electronic densities of states for the two models, 
as indicated. The identification of the peaks is discussed in 
the text. 

SSF's (or RDF's) can easily be achieved (cf. Figs. and 
^ , the interpretation of the small differences is not sim- 
ple, as they could be due to variations in the modes of 
preparation, details of the electronic potentials which 
could affect the structure without changing the con- 
nectivity, etc. The problem could however be resolved 
through measurements of the partial SSF's, as discussed 
above; unfortunately, this does not seem to be possible 
at present. 



V. ELECTRONIC PROPERTIES 

The electronic densities of states (DOS) for our two 
samples, CRN-A and CRN-B, are displayed in Fig. ||. 
These were obtained by averaging over the MD trajec- 
tories for a period of 3.5 ps at 300 K. As discussed in 
Refs. 25 and E6| the bands in the DOS can be roughly 



ascribed as follows: the lowest- lying band, labeled 'C in 
Fig. g, is As s-like while the next one, 'B', arises from 
Ga s and some As p states; the gap between this band 
and the following one is the "ionicity" gap. Just below 
the forbidden gap, band 'A' is composed of As p and Ga 
p states. In a crystal, the gap is direct and has a width 
of 1.55 eV. 

XPS measurements of amorphous GaAs have been in- 
terpreted as indicatipg that the material is essentially 
chemically ordered.cJ Probing the valence band from the 
gap down to about —15 eV, XPS reveals relatively lit- 
tle difference from the crystalline state except for the 
filling of the minimum between the first (A) and sec- 
ond (B) peak and a shift of 0.5 eV upwards of the third 



peak (C). However, upon comparing the DOS for our two 
networks (Fig. |6|), one sees an increased contribution in 
model CRN-A of the high-energy side of the low-lying 
As s band (C), almost forming an additional peak; this 
is due to wrong As- As bonds, as was also shown in Refs. 
psi and ^. The width of the gap at 300 K between this 
band and the following mixed band, —1.5 eV in the crys- 
tal, is reduced to about 0.9 eV in CRN-B and 0.5 eV in 
CRN-A. 

Another manifestation of the presence of wrong bonds 
is visible in the high-energy tail ofptiiG-|B band, which is 
much broader than in the crystal.lHal2a The added con- 
tributions at about -4.5 eV must be due to wrong bonds 
since the structural properties of the two models are very 
similar except for the number of such defects. This is in 
agreement with the discussion presented in Ref. Ell 

It is clear from the comparison between the two mod- 
els that although present, the effects of the existence of 
wrong bonds on the electronic structure of the material 
are much weaker than can hopefully be measured using 
techniques such as XPS. Thus, the observed similarities 
in the XPS spectra of a- and c-GaAs cannot be taken 
as evidepfe that the amorphous material is chemically 
ordered,c2l though our calculations do show that this is 
indeed the case. 

Disorder influences strongly the gap between valence 
and conduction bands. The value for this quantity has 
been rejiorted in the literature to lie in the range 0.61- 
1.45 eV.E!l Although the gap in the CRN-A sample is 
about 50% smaller than the one for CRN-B, both are 
substantially smaller than in the crystal. In particular, 
this value should depend significantly on the method of 
preparation of amorphous GaAs and is probably not a 
very good measure of the bond concentration.EEl 



VI. VIBRATIONAL PROPERTIES 

We have calculated the vibrational densities of states 
(VDOS) of our two models by Fourier transforming the 
velocity auto-correlation function averaged over 3.5 ps 
at 300 K; they are displayed in Fig. M. Experimentally, 
this quantity can be extracted from Raman spectroscopy 
measurements since amorphous networks do not have for- 
bidden symmetries so that all vibrational modes are al- 
lowed and measured. A multiple-order Raniap scattering 
study of a-GaAs has recently been reported.E^ The VDOS 
extracted from these data reveals two broad peaks, at 
about 2.3 and 8.3 THz, corresponding to the transverse 
acoustic and optic modes, respectively, with a wide, al- 
most featureless band in between.__At variance with the 
simulation results of Molteni et a/. Jl3 the TO peak is very 
much present experimentally, with a weight larger than 
that of the TA peak. 

From our results — Fig. ^ — we see that the VDOS for 
our two model networks are very similar: one difference 
is perhaps a slight shift of weight from the TA to the TO 
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FIG. 7. Partial and total vibrational densities of states 
for the two models, as indicated. Dashed lines are for Ga 
atoms, solid lines for As and dotted lines are the totals. The 
lower panel presents a comparison of CRN-A (dashed line) 
and CRN-B (solid line) with the crystal (dotted line). 

peak for CRN-B compared to CRN-A. The great sim- 
ilarity in the VDOS of our two models indicates that 
this quantity, also, is rather insensitive to the presence 
of wrong bonds; this would not be the case, however, 
if the two species differed appreciably in mass or elastic 
properties. In a recent Car-Parrinello simulation of a- 
InP, for instance, a high-energy peak in the partial PtE, 
VDOS has been identified as arising from wrong bonds;E2l 
because In is significantly heavier than P, the correspond- 
ing peak for In-In is lost in the "normal" continuum of 
states. 

Although the agreement between simulation and ex- 
periment is satifactory, a discrepancy remains regarding 
the relative amplitude of the two peaks. The weight 
of the TA peak is slightly larger than that of the TO 
peak in our simulation while the opposite is true exper- 
imentally. This difference indicates that "real" a-GaAs 
would have even less coordination defects that our model. 
Indeed, the TA peak is associated with bond-bending 
modes which are relatively insensitive to the local order 
as long as the stress is small, while the TO peak, depends 
critically on the e xiS it en ce of the tetrahedral symmetry 
around each atom;E2lE3 any coordination defect causes 
decrease of this peak. This is in fact clearly eviclent if we 
compare our results with those of Molteni et al.,t3 whose 
structure contains a much higher density of coordina- 
tion defects than ours and shows almost no TO peak. 
Based on this, therefore, and on the agreement of our 
calculated DOS with experiment, we conclude that real 
a-GaAs must be almost perfectly four-fold coordinated. 



VII. CONCLUDING REMARKS 

Using a newly-proposed, event-based, Monte-Carlo op- 
timization method, the activation-relaxation technique of 
Barkema and Mousseaup we have constructed a model of 
a-GaAs with a minimum of wrong bonds corresponding 
to an "infinite" Connell-Temkin model (CRN-B). This 
model is found to be energetically favorable over the tra- 
ditional Polk- type continuous random network (CRN-A). 
The CRN-B model of a-GaAs represents, to the best of 
our knowledge, the best realization to date of this mate- 
rial, with an almost perfect fourfold coordination, realis- 
tic bond-angle distribution, and almost no wrong bonds. 

In order to provide insight into the structure of a- 
GaAs, a detailed study of structural, electronic, and vi- 
brational properties of CRN-B has been presented, in- 
cluding a comparison with CRN-A. These results are 
in agreement with experiment and suggest that "real" 
a-GaAs forms a perfect CRN network, tetravalent and 
only weakly strained, with a minimum of wrong bonds 
(ideally none). Our analysis also shows, however, that 
wrong bonds are extremely difhcult to identify experi- 
mentally; in particular, indirect measurements (such as 
XPS and Raman scattering) cannot provide such infor- 
mation. Likewise, diffraction experiments that do not 
discriminate between the two chemical species are not 
sufficiently accurate to yield even approximate estimates 
of the proportion of wrong bonds. Our calculations indi- 
cate that only direct measurements of partial radial dis- 
tribution functions can provide experimental values for 
the density of wrong bonds. Such measurements, un- 
fortunately, do not seem to be possible at present. The 
results presented here thus provide a useful reference for 
further experimental and theoretical work. 
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